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Abstract 

Reduced models of neuronal activity such as Integrate-and-Fire models 
allow a description of neuronal dynamics in simple, intuitive terms and 
are easy to simulate numerically. We present a method to fit an Integrate- 
and-Fire-type model of neuronal activity, namely a modified version of the 
Spike Response Model, to a detailed Hodgkin-Huxley-type neuron model 
driven by stochastic spike arrival. In the Hogkin-Huxley model, spike 
arrival at the synapse is modeled by a change of synaptic conductance. For 
such conductance spike input, more than 70% of the postsynaptic action 
potentials can be predicted with the correct timing by the Integrate-and- 
Fire-type model. The modified Spike Response Model is based upon a 
linearized theory of conductance-driven Integrate-and-Fire neuron. 

Keywords: conductance injection - Integrate-and-Fire model - stochastic 
input - mapping techniques - predictive power. 

PACS: 87.10. +e - 87.19.La - 87.17.Nn - 87.17.Aa. 



1 Introduction 

The seminal work by Hodgkin and Huxley on mathematical description of 
action potential generation, has led to a whole series of models that describe 
in detail the dynamics of various ionic currents and the effect of the dendritic 
architecture on signal integration, see e.g. However, the precise 

description of neuronal activity involves a large number of variables, which often 
prevents a clear understanding of the underlying dynamics. Hence, a simplified 
description is desirable and has been subject of numerous studies (for a review, 
see ^112]). The most popular simplified models include the Integrate-and-Fire 
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model the FitzHugh-Nagumo model |5J|^] and the Morris-Lecar model |23) . 
However, it is not clear if such simplified models are sufficient to capture the 
essence of neuronal dynamics. Indeed, reduced models of neuronal activity, as 
opposed to detailed models of the Hodgkin-Huxley-type, are commonly thought 
to be too simple to account for the rich firing behavior of real neurons. 

Nevertheless, these highly simplified models have been shown to yield good 
predictions when compared to biological data [IT\. I n particular, several 
parameter estimation and/or optimization techniques have been proposed to 
map reduced models to real neurons ^] El E3 1201 HZ] • These techniques al- 
low to map electrophysiological data from biological neurons (intracellular or 
extracellular recordings) onto simplified models of the fntegrate-and-Fire-type. 
Moreover, such a mapping could be the starting point of a systematic classi- 
fication of cortical neurons in terms of simplified dynamics. Finally, reduced 
models have been used extensively and successfully to model and study analyt- 
ically the behaviour of cortical- like networks 0J^2- I n other words, analytical 
and computational tools are available to go from a simple description of single 
cells dynamics to network dynamics. 

In this paper, we review and extend a numerical technique that allows a 
systematic mapping of a class of Integrate-and-Fire neuron models, namely the 
Spike Response Model ^JE]> to intracellular recordings of neuronal activity 
|T5J. While this technique ensures reliable predictions and good generalization 
when the target neuron is driven by a randomly fluctuating current, it was 
shown in a previous study |15j that the generalization power is poor when the 
target neuron is driven by randomly fluctuating synaptic conductances. Hence, 
the purpose of this paper is to propose extensions of our previous work which ad- 
dresses this latter problem. Using recent theoretical results [20 CHI > we propose 
a generalization of the Spike Response Model so as to make model parameters 
input-dependent. This improved version of the SRM is shown to be very effi- 
cient and robust at predicting the spike train of a detailed Hodgkin- Huxley type 
neuron model. 

2 Materials and Methods 
2.1 The Spike Response Model 

We consider a neuron stimulated by stochastic presynaptic spike arrival. The 
state of the neuron is characterized by a single variable u, the membrane voltage 
of the cell at the soma. Let us suppose that the neuron has fired its last spike 
at time t. At time t > t, the membrane potential of the cell is described as: 

u(t) = V (t -i)+EE e+ (*- *f) +EE r <*- *?)• c 1 ) 

i£E f jei f 

The last two terms account for the drive by presynaptic neurons i (respectively 
j) of the excitatory E (respectively inhibitory I) population, tj and t- denote 
the firing time of presynaptic neurons. The e-functions describe the excitatory 
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(+) and inhibitory (— ) postsynaptic potentials (EPSPs and IPSPs). rj(t — i) 
describes the form of the spike itself as well as the after-hyper-polarization 
potential, if present. A spike is elicited if the following threshold condition is 
satisfied: 

if u(t) > and — > then, t = t. (2) 

Note that spiking occurs only if the membrane voltage crosses the threshold 
fl from below. The threshold itself can be taken either as a constant or as 
time-dependent. In this paper, we use a dynamic threshold with the form: 

»/. _ a _ f +00 if t - 1 < lref , . 

[ ' \ ?9o + ??iCxp(-(<-f)/r 1? ) else ' [ ' 

where 7 re / is a fixed absolute refractory period so as to exclude continuous 
firing. #Qj $i an( A tq are parameters that will be chosen to yield the best fit to 
a target spike train (see section |2~21 below) . This version of the Spike Response 
Model (the one that we use in this paper) is a simplified version of the full Spike 
Response Model and has been termed SRMo (see for further details). We 
will use this acronym to refer to this model. 
Eq. (JJJ can be restated in the form: 

r+oo r+oo 

u(t)=7](t-i)+ e+(s)Q+(t- s)ds+ e-(s)Q-(t-s)ds (4) 



with: 

Q + (*) = EE^-*/); ^(*) = EE^-*J)- ( 5 ) 

/ jei f 

For numerical implementation, we will use a discrete version of Q's: 

r t+At 



Q + -{s)ds (6) 

i.e. Qf is the spike count in a time bin of duration At in the excitatory presy- 
naptic population and analogously for . The activity (population rate) A + 
of the excitatory presynaptic population is defined as follows: 

where N + is the size of the presynaptic population (with a corresponding defi- 
nition for the inhibitory population). 

2.2 Mapping the SRM to voltage traces 

The mapping technique that we propose has been discussed in detail elsewhere 
[T5] and we refer interested readers to this reference. Here, we describe the 
essentials of the technique without going into details. 



3 



To realize the mapping of the SRMo to the target neuron, we proceed in two 
steps. First, we extract the functions characterizing the model (EPSP e + , IPSP 
e~ and spike shape 77) and second, we choose the parameters of the dynamic 
threshold (1?) and optimize them in terms of quality of predictions. To do so, 
we assume that we have at our disposal voltage traces of the target neuron 
as well as firing times of presynaptic neurons. We also assume that the input 
characteristics are kept constant during the recording of the dataset used for 
the mapping procedure. We start by extraction of the spike shape 77. The shape 
of spikes is usually highly stereotyped and presents only little variability. We 
therefore select one spike train from the dataset and align all spikes relatively to 
some arbitrarily chosen initiation point. The mean trajectory of the spikes yields 
77. Detection and alignment of spikes is realized using a threshold condition on 
the first derivative of the membrane voltage. Once we are done with 77, we 
extract the shape e + of an EPSP and the shape e~ of an IPSP. If we limit 
ourselves to the interval between two consecutive spikes of the same spike train 
tk and ifc+i, we can rewrite Eq. with the notation introduced in Eqs. 
and © as follows (for ij~ < t < tfe+i): 



u(t)-r](t-tk)= e + (s)Q+(t-s)ds+ e-(s)Q-(t-s)ds. (8) 



It is then possible to find the optimal e-functions using the Wiener-Hopf optimal 
filtering technique 20, 32] . We fit the resulting EPSP e + with a suitable function, 
typically a difference of exponentials: 



with a corresponding definition for e - . K e , r r and are free parameters. The 
final step is to choose and optimize the threshold. The absolute refractory pe- 
riod "/ re f is set to 2 ms. The other parameters of Eq. i.e. t?oi $1 an d t#, are 
fitted in order to optimize the coincidence factor T (see below) on a given spike 
train. In order to ensure a good generalization, we optimize the threshold with 
a large dataset generated with different input characteristics. The numerical 
optimization algorithm that we use is the downhill simplex method |25|. Obvi- 
ously, the SRMo can only predict neuronal activity of the specific neuron it has 
been mapped to. 

2.3 Target neuron model 

Instead of real data from experiments, we use as a reference or "target" a 
Hodgkin- Huxley- type model of a fast-spiking interneuron [5]. It contains stan- 
dard Na + and K + spike-generating channels and one extra Kvl.3 K + -channcl 
derived from "n"-type currents measured in human T-lymphocytes. This latter 
current produces a subtle form of adaptation. See |S| for further details. 





e+(s) = K+ (exp(-s/r+) - exp(-s/rj 



))■ 



(9) 
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synapse 


E (mV) 


T syn (ms) 


D (mS/cm 2 ) 


exc. (+) 





2.45 


0.073 


inn. (-) 


-80 


6.11 


0.04 



Table 1: Parameters of excitatory and inhibitory synapses (adapted from Des- 
texhe and Pare [H]). Parameter D has been adjusted so as to yield an amplitude 
of postsynaptic potential of the order of 1 mV. 

2.4 Input scenario 

The Hodgkin-Huxley-type model of a fast-spiking interneuron is driven by ran- 
dom synaptic conductances generated by massive stochastic presynaptic spike 
arrival. The total driving current is given by |29| : 

W*) = 9 + (t) (u(t) - E+) + g-(t) (u(t) - £_), (10) 

where g+,~ are the excitatory (+) and inhibitory (— ) conductances and 
are the corresponding reversal potential. Synaptic conductances are generated 
by slightly correlated homogeneous Poisson spike trains. The excitatory and 
inhibitory populations contain respectively N + = 8000 and N~ = 2000 neurons 
of which only a fraction N + /N + (respectively N~ /N~) are independent. N + '~ 
depend on the presynaptic discharge frequency, on the correlation coefficient 
and on the size of the presynaptic population. In the following, we assume that 
the presynaptic neurons discharge at frequencies v + and v~ and the correlation 
coefficient is held at a constant value c = 0.002 (see Appendix for further 
details about how the presynaptic spike trains are generated and, in particu- 
lar, see Eqs. (|31(l and i|32|) for the relation between N + '~ and the correlation 
coefficient c). 

The dynamics of each excitatory synapse i is described by a conductance 
variable with: 

n+ e+oa 

t+ d P+ = -P+ + D+£<*(t-tf) = — / e- s ^Q+(t~s)d S . (11) 

p Tsyn J 

The value of P^~ is increased by an amount D + for each presynaptic spike acti- 
vating the synapse at time ti. It then decays back to zero with a time constant 
r+ n . The total excitatory conductance is the sum of conductance variables 
over all excitatory synapses. Thus, the total excitatory conductance is: 

<4a(*) = E^ + (*)> (12) 

with a corresponding definition for the total inhibitory conductance g^y n (t) = 
y^jef Numerical values used in this paper are summarized in Tabled 
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2.5 Assessing the quality of predictions of the reduced 
model 



In order to evaluate quantitatively the predictions of our reduced model, we use 
the coincidence factor T between two spike trains ^7] defined as: 



where -/V ta rget is the number of spikes in the target spike train Starget, -^SRM is 
the number of spikes in the spike train 5s rm that is predicted by our reduced 
model, iVcoinc is the number of coincidences with precision A between the two 
spike trains, and (iV CO i nc ) = 2^AiVtargct is the expected number of coincidences 
generated by a homogeneous Poisson process with the same rate v as the spike 
train 5srm- The factor TV = 1 — 1v A normalizes T to a maximum value of one 
which is reached if and only if the spike train of the reduced model reproduces 
exactly that of the target neuron. A homogeneous Poisson process with the 
same frequency as the reduced model model would yield, on average, T = 0. 
We compute the coincidence factor F by comparing the two complete spike 
trains, i.e., the spike train Starget generated by the target neuron and the train 
■S'srm predicted by the SRM. This is different to the approach of Kistler and 
colleagues |17| where T was used to predict the next spike in a spike train, under 
the assumption that past action potentials were correctly reconstructed. 

2.6 Linearized theory of a conductance-driven Integrate- 
and-Fire model 

The results of section 3.1 show that a SRMo with fixed time course of EPSPs 
e + and IPSPs e~ has a rather limited range of validity. The naive solution 
would therefore be to use a different set of EPSPs and IPSPs for each set of 
discharge frequencies However, this is not a very practical solution 

if we have to compute PSPs e + and e _ for each possible combination of input 
frequencies with the method indicated above. Instead, we aim at finding a sim- 
ple parameterization of the EPSPs and IPSPs that would allow to interpolate 
between and generalize beyond the specific inputs used for the numerical deriva- 
tion of e + and e~ . To do so, we use a linearized theory of conductance-driven 
Integrate-and-Fire models that allows to write down an analytical expression 
for the EPSPs e + and IPSPs e~ in function of the discharge frequencies v + and 
v~ . Thus, the extended Spike Response Model SRM C that we derive in this 
paper is directly related to the subthreshold dynamics of a conductance-driven 
Integrate-and-Fire neuron model (CIF) [51 1281 l3*U| . 

The subthreshold membrane voltage u of a CIF neuron model is given by 
the following differential equation: 



r 




(13) 



2 (-^target + N SB , U ) AT 



r d 
C —u 

dt 



g L (u - E L ) ~ g+(t) (u - E+) - g~(t) [u - E~), 



(14) 



G 



with C the membrane capacitance, <?l the leak conductance (with a reversal po- 
tential El) and g + (respectively g~) the instantaneous excitatory (respectively 
inhibitory) conductance. E + and E~ are the reversal potentials of the exci- 
tatory and inhibitory synapses. It is straightforward to show that the average 
membrane voltage predicted by this model (in absence of spiking mechanism) 
is given by: 

gLE L +gt n E + + g- yn E^ 
Mcif = , _+ , — ■ (15) 

9L + Ssyn + 3syn 

The bars denote time averaging in this case. The EPSPs e + and IPSPs e~ of 
the SRMo that would correspond to the CIF model can then be calculated as 
the linear response around the average membrane voltage ^cif- We find: 

D+ t cS r+ n (/xcif - E+) / + _ eff - 



e+ (s) = hyl lvr 7 - e- s ' T ^ - e- s l T ^ . (16) 

T e S — T syn ^ ' 

with a corresponding definition for e~(s). We compare the results of Eqs. I|15fl 
and (|16|l with results extracted from simulations of the target Hodgkin-Huxley- 
type neuron model. We find that, outside the regime of EPSP-amplification, 
both the average membrane voltage (data not shown) and the shape of PSPs e + 
and e~ are well predicted by this simple linearized theory (see Figure|2K and B 
in Results section). To illustrate this, we fit the PSPs plotted in Figure^V with 
Eq. (|16|) with r e ff as a free parameter. We then compare this parameter with 
r e ff as predicted by the Eq. (|19(l . Figure |3] shows that the linearized theory is 
in very good agreement with the results of simulations of the target Hodgkin- 
Huxley-type neuron model except in the area where EPSP-amplification takes 
place. 

Equation l|16|l gives an analytical expression for the PSPs e + and e~. The 
PSPs allow us to reproduce the fluctuations of the membrane voltage. How- 
ever, we also need to account for the constant voltage bias which appears when 
discharge frequencies are elevated. Therefore, we redefine the first term of the 
SRMo, namely the function r/ as: 

r)(t-t)^r)(t-t) + (ticiF-E L ). (17) 

This simple procedure ensures that our model produces the correct average 
membrane voltage. The model could be further improved by using a time- 
dependent leak conductance gL {t — i) |15| , but we will not do so. 

In order to relate the linearized theory of conductance-driven Integrate-and- 
Fire model of section 3.2 to the numerical PSPs e + and e~ extracted in sub- 
section 3.1, we need to estimate the parameters appearing in Eq. I|16|) from 
the data. We will assume that the size of the presynaptic populations (iV + 
and N~) as well as the average discharge frequencies in these populations {v + 
and v~) are known. We will also assume "standard" reversal potentials at the 
synapses, i.e. E + = OmV and E~ = — 80 mV. C is taken to be constant at a 
value of 1/xF/cm 2 . The reversal potential of the leak current can be crudely 
approximated by the resting potential of the neuron El ~ u rcst . 
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Thus, the parameters we need to estimate are the synaptic time constants 
T s + n and r s ~ n , the effective membrane time constant r c g, as weil as the mean 
conductances gL, g + and g~. The time constants r+ n (respectiveiy T s ~ n ) can be 
estimated by fitting the EPSP e + (respectiveiy the IPSP e~) extracted from a 
dataset where EPSP-amplihcation doesn't take place, i.e. a dataset with strong 
inhibition and weak excitation. Once T+ n and r S y n are known, one can find 
the other parameters by comparing EPSPs e + extracted for several different 
sets of input characteristics. Fitting the EPSPs e + with formula (|16|1 yields an 
estimate of r c g. If we collect three different values of r c g for three different sets 
of input characteristics {v + and and using the definition of r e ff (see Eq. 
I|19l) in Results section), we find a set of three equations with three unknowns, 
namely g L , g+ n and <7" n . g+ n is given by: 

g+ n = N+v+T+ n D+ (18) 

with a corresponding definition for <?^ n . Solving this set of equations thus yields 
an estimate of gL, D + and D~ . Therefore, we now have a simple model that 
allows to compute PSPs e + (s) and e~(s) in a straightforward fashion for any 
given set of input discharge frequencies v + and v~ . We will refer to this model 
as SRM C with a lower-case "c" for conductance-based. 

3 Results 

In the first subsection, we discuss the intrinsic limitations of the simple Spike 
Response Model (SRM ). In particular, we show that the shape of EPSPs 
and IPSPs derived by our method changes systematically as a function of the 
input parameters. These findings are then compared with a linearized theory of 
conductance-driven Integrate-and-Fire models (CIF). This comparison allows us 
to determine the parameters of a new conductance-based Spike Response Model 
(SRM C ) which is tested over a broad range of different inputs in subsection 3.2. 

3.1 Limitations of a classic SRM 

We map the target neuron model to the SRMo using the technique described in 
Materials and Methods. Let us recall that the SRMo is characterized mainly by 
the spike shape r\ as well as the EPSP e + and the IPSP e~. Fig. shows the 
EPSPs and IPSPs e + '~ extracted for different input discharge frequencies v + 
and v~ and Fig. Qj3 shows the spike shape rj. The shape of the spike does not 
depend on the characteristics of the input scenario. One immediately remarks 
that both the characteristic time scales and the amplitude of the EPSP/IPSP 
do change in function of the input discharge frequencies. In fact, the numerical 
technique extracts the best linear filters so as to reproduce the membrane voltage 
trace of a given sample spike train. The corresponding EPSPs and IPSPs are 
then optimal for the specific set of input discharge frequencies used for parameter 
extraction but there is no reason why they should be optimal for other sets of 
inputs with different characteristics. Indeed, there are reasons why they should 
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be different depending on input discharge frequencies. Let's quickly review these 
reasons [H]. 

First, two different sets of input discharge frequencies are likely to produce 
two different values of average membrane voltage. In turn, this means increased 
or decreased average driving forces for the synapses as the corresponding cur- 
rent includes a multiplicative term of the form (u(t) — E + '~) (see Eq. 
This should affect the amplitude of the EPSPs and IPSPs. Second, two differ- 
ent sets of input discharge frequencies are likely to produce two different total 
conductances. This affects the effective membrane time constant of the neuron. 
If we neglect the effect of somatic AP-generating ion channels, the effective time 
constant in the subthreshold regime can be written: 

Toff - + , (19) 

9L + glyn + ffsyn 

with C the membrane capacitance, gh the conductance of the leak current and 
g^ yn and g~ yn the total average excitatory and inhibitory synaptic conductances. 
Finally, when the excitatory input discharge frequency is very high, the target 
neuron runs in a highly non-linear regime due to activation of Na + -channels. 
This effect is known as EPSP-amplification ^U] and illustrated for our target 
neuron in Fig. \I]p and D. While EPSP-amplification is usually not observed 
in fast-spiking interneurons, it is present in our target model of a fast-spiking 
interneuron since Na + -channels are the only channels opened at depolarized 
states close to threshold (see Materials and Methods section). All these three 
effects (average driving force, effective membrane time constant and EPSP- 
amplification) combine with each other and lead to the pattern of EPSPs and 
IPSPs shown in Fig. [TJA and C. The EPSPs and IPSPs are shorter when total 
conductance is increased (shortening of effective membrane time constant) ex- 
cept when this increase is mainly due to the excitatory conductance which then 
leads to the reverse effect because of EPSP-amplification (see panel C). The 
amplitude of EPSPs decreases when excitation is increased and increases when 
excitation is decreased (average driving force) with corresponding effects for IP- 
SPs and modifications of the inhibition level. Interestingly, one remarks that 
when inhibition is very strong together with weak excitation, the extracted EP- 
SPs e + and IPSPs e~ follow exactly the dynamics of the corresponding synapses, 
i.e. an almost instantaneous rise followed by a decay with the same time con- 
stant as the time constant of the synaptic conductance (see Fig. at bottom 
right). The dynamics of EPSPs and IPSPs can be approximated by: 

e + -(s) oc (e- s/ <* - e~ s/T ^ . (20) 

Therefore, when drive is very strong, r e ff ~ from Eq. Ijl9[) and e +, ~ is reduced 
to e~ s / rs y" . Note that this effect does appear only with weak excitatory stim- 
ulation. When excitatory drive is strong, EPSP-amplification compensates the 
reduction of T e s and Eq. (|19fl does not hold anymore. 

Given the change of the time course of EPSPs and IPSPs as a function of 
the input, it is clear that correct predictions of the subthreshold fluctuations of 
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Figure 1: EPSPs, IPSPs and spike shape as extracted by our numerical method. 
A. EPSPs and IPSPs (respectively positive and negative solid lines) for differ- 
ent presynaptic input discharge frequencies v + = 0.9, 0.6, 0.3 Hz (from top to 
bottom) and v~ — 1, 2, 3, • • • , 12 Hz (from left to right). Horizontal axis is in 
ms and vertical axis is in mV. The light grey area shows the dynamics of the 
normalized EPSCs and IPSCs (vertical axis in arbitrary units). B. The spike 
shape r\ C. Effect of Na + -channels on the shape of kernels. This panel repeats 
the results of panel A with discharge frequencies v~ = 5 Hz and v + = 0.9 Hz 
(light grey area and solid lines). The dotted lines correspond to EPSPs and 
IPSPs extracted while blocking Na + -channels. Resulting PSPs are significantly 
shorter. D. EPSP-amplification in the target model. An exponentially decaying 
excitatory conductance is injected in the target neuron (bottom; vertical bar is 
0.005 mS/cm 2 ). The response of the target neuron (top) when receiving n-fold 
times the excitatory conductance for n = 0,4, 8, 12, 16, 20 (thin solid lines). For 
n = 20, one observes huge EPSP-amplification in normal conditions (thick solid 
line) but none when Na + -channels are blocked (thick dotted line; note that in 
this latter case, the resting state is slightly changed) and a spike is finally elicited 
for n — 21. Horizontal bar is 10 ms and vertical bar is 10 mV. 
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the membrane voltage by linear summation of EPSPs and IPSPs with a fixed 
time course are only possible in the input regime for which the reduced model 
was constructed (recall Eq. QJ). 

In order to make the EPSPs and IPSPs input dependent, we parameterize e + 
and e~ using the linearized theory of a conductance-driven Integrate-and-Fire 
model 5, 28, 3Uj; see Materials and Methods for details. 

We find that outside the regime of EPSP-amplification the numerically de- 
rived postsynaptic potentials e + and e~ are well fitted by the theory (see Figure 
|2J). The regime of EPSP-amplification can easily be identified by comparing the 
effective membrane time constant T c g predicted by the theory with that derived 
from the numerically optimized PSPs e + and (see Figure yj. All parameters 
of the linearized theory, in particular the synaptic time constants r+ n , r s ~ n and 
the mean conductances g + and g~ can hence be estimated from the data by 
using three sets of inputs that do not lead to EPSP-amplification (see Materials 
and Methods). The resulting model is a conductance-driven Spike Response 
Model (SRM C ) which we now test on new set of inputs independent from the 
one used for parameter optimization. 
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Figure 2: The EPSPs and IPSPs as predicted by the linearized theory (solid 
lines; see Eq. (|lCi|> ^ are compared to the EPSPs and IPSPs extracted by the 
method proposed in Materials and Methods (symbols). A. The inhibitory 
discharge frequency v~ = 6 Hz and the excitatory discharge frequency v + = 
0.9, 0.6, 0.3 (from top to bottom). B. Same as in A except that v~ = 10 Hz. 
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Figure 3: The effective membrane time constant as predicted by a C1F (soiid 
line) is compared to the corresponding parameter extracted from simulations of 
the target Hodgkin-Huxley-type neuron model (dotted line with squares). Sec 
text for further details. A. v+ = 0.3 Hz. B. v+ = 0.6 Hz. C. v+ = 0.9 Hz. 
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3.2 Predicting spike by spike 

A conductance-driven Spike Response Model (SRM C ) has been estimated from 
the numerical voltage traces using the procedure described in the Materials 
and Methods section. We now test the predictive power of the SRM C . We are 
interested at reproducing the exact timing of the spikes of the target neuron. As 
we have not dealt yet with the threshold, the first step is to optimize the three 
free parameters of the threshold, namely i?o, #1 and t# (see Eq. ©, 7 re / is set 
at a constant value of 2 ms). To do so, we use a very long spike train with input 
characteristics that include the discharge frequencies where EPSP-amplification 
doesn't take place (see Figure EJ. We then proceed as indicated in Materials 
and Methods. The resulting numerical values are summarized in Table [5] 



parameter 


mean 


SD 


tfo (mV) 


-38.437 


0.002 


01 (mV) 


564.0 


0.7 


Ttf (ms) 


0.91 


0.03 



Table 2: Best fit parameters for the threshold (see text for further details). 
Mean and standard deviation (SD) are computed from four optimizations with 
different initial conditions. The small standard deviation shows that the four 
optimization runs all converge to the same minimum. 

Figure 0] shows sample results of the predictive power of the SRM C in two 
distinct regimes, low-drive (low presynaptic discharge frequencies, 1) and high- 
drive (high presynaptic discharge frequencies, 2). In both cases the predicted 
membrane voltage is reasonably close to the membrane voltage of the target 
neuron. Note also that all the spikes in plotted samples are reproduced with 
the exact timing (panel B). For these two cases, we find V = 0.76 (low-drive) 
and r = 0.67 (high-drive). The coincidence factor T takes a value of 1 if 100% of 
spikes coincide and is normalized to if coincidences are random (see Materials 
and Methods section). 

To test the performances in a more systematic way, we quantify the pre- 
dictions of the SRM C in terms of the timing of the spikes (coincidence factor 
F) and in terms of output frequencies of the model (f ou t) over a broad range 
of input characteristics (v + and v~). Figure shows the performances of the 
SRM C for such a systematic procedure. We observe that the SRM C yields good 
performances (r > 0.7 and output frequency v avAi predicted in the correct range) 
over a broad range of input discharge frequencies, the only notable exception 
is when inhibitory discharge frequency is very high. In this regime, the output 
frequency v out is close to zero so that the number of spikes in a spike train is low. 
Therefore, our coincidence factor T is not well suited for this case. However, 
we note that even if T is low, the output frequency is predicted in the correct 
range. Furthermore, the subthreshold fluctuations of the membrane voltage are 
well reproduced (not shown). 
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Figure 4: Comparative results for low (1) and high (2) input discharge frequen- 
cies. A. The activity in the excitatory (top) and inhibitory (bottom) presynap- 
tic populations. The arrows indicate the starting point of the segment plotted 
in B. B. Corresponding membrane voltage of the target neuron (solid line) is 
compared to the membrane voltage as predicted by the SRM C (dashed line). In 
both cases (B.l and B.2) the membrane voltage as predicted by the SRM C gives 
a fair approximation of the membrane voltage of the target neuron. 



Interestingly, the SRM C yields good performances where it is supposed to 
do so but also where EPSP-amplification takes place (see Figure 0. In the 
case of our target neuron, EPSP-amplification occurs when both excitatory and 
inhibitory populations discharge at rather low frequencies. A fixed coefficient of 
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correlation c then imposes co-activation of large subpopulations of synapses in 
this regime (see Eq. H310 ). The target neuron therefore spends most of the time 
close from the resting potential and occasionally receives numerous coincident 
presynaptic spikes that bring it in the area of EPSP-amplification. Most of 
these large excursions then lead to a spike while only a few of them do not. In 
this situation, one can easily predict spikes at a correct timing with a threshold 
model like the SRM C . 
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Figure 5: Coincidence factor T (A; mean ± SD) and output frequency v OVL t (B) 
plotted versus the inhibitory input discharge frequency v~ for various values of 
the excitatory input discharge frequency v + (v + = 0.3 Hz solid line and squares; 
v + = 0.6 Hz dashed line and diamonds; v + = 0.9 Hz dotted line and circles). In 
the panel B, the output frequency of the SRM C (symbols) is compared to the 
output frequency of the target neuron (dots and line). 

In summary, our modified version of the SRMo, the so-called SRM C , performs 
well over a broad range of input characteristics v + and v~ and predicts not only 
the output frequency but also most of the spikes with the correct timing. It is 
useful to keep in mind that we used a calculation based on a conductance-driven 
Integrate-and-Fire model to evaluate the subthreshold membrane voltage fluc- 
tuations and the average membrane voltage. This latter results do not include 
the potential effect of ionic channels activated in the subthreshold regime. The 
SRM C should therefore fail in regimes where such ionic channels are activated. 
However, we found here that even though the underlying theoretical framework 
is on the edges of its validity domain, the SRM C still performs reasonably well 
in the regime where EPSP-amplification takes place and thus as long as there 
are large enough fluctuations of the subthreshold membrane voltage. 
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4 Discussion 



Mapping real neurons to simplified neuronal models has benefited from many 
theoretical developments in recent years and applied to both in-vitro and in-vivo 
recordings |161 177| . However, most of the techniques have been developed for a 
current injection scenario [51 El 1261 127| . On the experimental side, conductance 
injection is increasingly used instead of current injection and is thought to be 
closer to in-vivo conditions (see 7 for a review). 

We had previously reported a mapping technique |15j based on standard 
signal processing tools which allows a systematic mapping of a simplified neuron 
model, the Spike Response Model to intracellular recordings. It has 

been shown to yield very good results in the case of current injection for model 
neurons |15j and with real data (unpublished observations). However, while the 
reduced model built in this way generalizes its predictions over a broad range 
of different input characteristics for the current injection case, it performs very 
poorly in the conductance injection case for reasons reviewed above. 

In this paper, we have shown that a simple modification of the classic SRMo 
solves this problem. In its new formulation, the model is able to predict very 
reliably many aspects of neuronal activity, such as timing of the spikes, mem- 
brane voltage and mean output rate. The global performances are improved 
and moreover, the model can generalize predictions extremely well. The new 
model SRM C is directly related to conductance-driven Integrate-and-Fire neu- 
rons jSJ |551 1301 • Our technique can be applied to extract simple neuron models 
from experimental intracellular recordings under conductance injection. 

A Generation of presynaptic spike trains 

In this appendix, we detail the method used in simulations to generate slightly 
correlated spike trains and we also derive some useful analytical results. The 
method follows |Hj but there exists other ways to generate correlated spike trains 
(see ^B] for instance). 

Presynaptic spike trains are described by random homogeneous Poisson pro- 
cesses. At each time step, N independent random variables are generated and 
distributed among the N > N presynaptic neurons to generate slightly cor- 
related spike trains |6J. In this specific case, we can derive the probability 
distribution function (PDF) of the variable Q (see Eq. © in Materials and 
Methods section for a definition of variable Q). Here, all the calculations rely 
on a discrete time scale with bins of width At. Let us consider that elements 
of a vector V of length N are distributed at random and receive either a value 
of 1 with a constant probability p or with probability (1 — p). In this spe- 
cific case, p = v At with v the average discharge frequency in the presynaptic 
population and At the size of the time steps used in the simulation. The total 
number K of elements receiving a value of 1 in V is therefore distributed ac- 
cording to a binomial distribution P(K = k) = B(k; N ,p). In a second vector 
V of length N, elements receive or 1 according to a parent element chosen at 
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random in vector V. The probability of receiving a value of 1 is then p — K/N. 
The total number of elements Q of V receiving a value of 1 is then given by 
P(Q = q\K) = B(q; N,p). The average distribution of variable Q is then: 



P(Q = q)=J2 P (Q = l\K)P(K = k )- (21) 



fe=0 

Some algebra yields: 

N 

P(Q = q) = Cf N~ N c k P k (1 - p f~ k k q (N - k) N - q . (22) 

fc=0 

In the following, we will need to know the first two moments E[Q] and Var[Q] 
of this distribution. Using the definition of p (see above), we find that E\p] = p 
and thus: 

E[Q] = Np. (23) 
To calculate the variance of Q, we use the fact that: 

E[Q 2 } = Var[Q]+E[Q} 2 

= Np(l-p)+N 2 f (24) 

Using the fact that E[p] = p and that E[p 2 } — N~ 2 Var[K] + p 2 , we find that: 

Var[Q] = Np (1 — p) (l + N/N) (25) 

While variables N and p have direct biological interpretations (N is the size 
of the presynaptic population and p is related to the discharge frequency in that 
population), TV is a rather abstract quantity which is linked to the correlations 
in the activity of the presynaptic population. For practical use, it would be 
useful to compute the correlation coefficient between two spike trains in terms 
of these variables. The correlation coefficient between two sequences of numbers 
nii and rij is defined by: 



E 4 ™f E 



(26) 



If we choose at random two elements from the vector V, the probability that 
both fire together is given by the hypergeometric distribution H (N, q, 2) so that: 

P(Y = m = ^ = 0A. y (27) 

On the other hand, if we choose only one element from the vector V, the prob- 
ability that it fires is given by H(N, q, 1): 

P(Z = l\Q)=q/N (28) 
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Finally, the average correlation coefficient between two randomly chosen spike 
trains of the presynaptic population is given by: 

c = E^oP(X = 2\Q)P(Q = q ) 
Y% =0 P(Z = l\Q)P(Q = q>)' 

Using Eqs. J23 and (J2EJ in Eq. we find: 

1 [Q 2 ] - E[Q] 
N-l E[Q] 



(30) 



Note that this latter result is general when considering homogeneous Poisson 
spike trains with the same rate p and doesn't depend on the specific way spike 
trains are generated. In our case, some algebra yields: 

, N(l-p) 

C=P+ N{N^T) (31) 

and thus: 

_ N{l-p) 

N -(N-i)(c- P y (32) 
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